theta = linspace(0, 2*pi, 1000);
z = (exp(1i*theta) - 1)./exp(1i*theta);
plot(real(z), imag(z)); 

hold on
z =  (3*exp(2 * 1i*theta) - 4*exp(1i*theta)+1)./(2*exp(2*1i*theta));
plot(real(z), imag(z)); 

hold on
z = (11*exp(3 * 1i*theta) - 18*exp(2 * 1i*theta)+9*exp(1i*theta)-2)./(6*exp(3*1i*theta));
plot(real(z), imag(z)); 

hold on
z = (25*exp(4 * 1i*theta) - 48*exp(3 * 1i*theta) +36*exp(2 * 1i*theta)-16*exp(1i*theta)+3)./(12*exp(4*1i*theta));
plot(real(z), imag(z)); 